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ABSTRACT 

We measure the rotation of the Milky Way stellar halo on two samples of Blue Hori- 
zontal Branch (BHB) field halo stars from the Sloan Digital Sky Survey (SDSS) with 
four different methods. The two samples comprise 1582 and 2563 stars respectively 
and reach out to ~ 50 kpc in galactocentric distance. Two of the methods to mea- 
sure rotation rely exclusively on line-of-sight velocities, namely the popular double 
power-law model and a direct estimate of the de-projected l.o.s. velocity. The other 
tw o techniques use the full 3D motions : the radial velocity based rotation estimator 
of ISchonrich. Binnev fc Asplundl ()2012l ) and a simple 3D azimuthal velocity mean. 
In this context we a) critique the popular model and b) assess the reliability of the 
estimators. All four methods agree on a weakly prograde or non-rotating halo. Fur- 
ther, we observe no duality in the rotation of sub-samples with different metallicities 
or at differen t radii . We trace the rotation gradient across metallicity measured by 
iDeason et al.l ()2011[ ) on a similar sample of BHB stars back to the inclusion of regions 
in the apparent magnitude-surface gravity plane kno wn to be contaminated. In the 
spectroscopically selected sample of IXue et al.l (|201l[ ). we flag ~ 500 hot metal-poor 
stars for their peculiar kinematics w.r.t. to both their cooler metal-poor counter-parts 
and to the metal-rich stars in the same sample. They show a seemingly retrograde 
behaviour in line-of-sight velocities, which is not confirmed by the 3D estimators. 
Their anomalous vertical motion hints at either a pipeline problem or a stream-like 
component rather than a smooth retrograde population. 
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1 INTRODUCTION 

The stellar halo represents an observable tracer of the 
putative dark matter halo which is thought to gov - 
ern our cosmological evolution (e.g. IWhite fc Reeal201(J ). 
The present observational evidence favours a hierarchi- 
cal formation scenario, where the stellar ha lo was formed 
through accretion of protogalact ic fragments l|Searle fc ZinnI 

/ after the detec- 
tion of abundant structures with SDSS (jBell et al.l 120081 : 



throu gh 

[lizi;! 



tog 



Bullock fc Johnston! 120051 '). especiall y after the detec- 



iBelokurov et al.l 120061 ). Simulations of stellar halos support 
the picture of a hal o mos tly made of accr eted structures (e.g. 
ICooper et al.ll201(]| '). but iFreemanI (|201ll ') remarks that a for- 
mati on process via dissipation cannot be c ompletely ruled 
out l|Eggenlll962l : ISamland fc Gerhardll20ol '). 

A transition from a flattened inner to a more spherical 
outer halo has sometimes been traced back to halo forma- 
tion in different stages and via different physical processes 
l|Hartwicklll987l : [Preston et allllQQl. '). Since it has been ar- 
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gued that dynamical friction is more eflicient fo r prograde 
infall (|Quinn fc GoodmanllToS^ : iBvrd et al.ll 19861 ) there is a 
popular expectation for the outer halo to have a d ifferent ro- 
tational signature from the inner halo IIMurante et al. 2010 ). 
and in particular to be counter-rotating (|Norriall989l ). 

The kinematics of the stellar halo could serve as a 
probe of formation history, but it is still highly contro- 
versial: on the one hand there are claim s for a globally 
counter-rotat i ng ha l o (e.g . iMaiewskil Il992). More recently 
ICaroUo et al.l (|2007l . [2010|) using solar neighbourhood main 
sequence stars claimed a inner pro-rot ating halo opposed to 
an ou ter counter-rotating halo, which iKinman et al.l (|2007l . 
[2013) conflrmed with a small sample of RR Lyrae stars and 
blue horizontal branch stars localized near the Galactic an- 
ticentre direction. With a larger sample of blu e horizon- 
tal branch field halo stars iDeason et al.l (|201ll ) found in- 
stead a metal-rich pro-rotating halo opposed to a metal-poor 
counter-rotating halo with such a dichotomy being coh erent 
at all radii . On t he other hand (iChiba fc BeersI [2000|) and 
ISirko et al.l (|200J) argue for a non-rotating halo, but do not 
distinguish for different metallicity. 
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The uncertainty in this debate arises from several 
causes: the challenge of a clean selection of halo stars, the 
question if such samples give a fair representation of the 
halo (e.g. still debated for globular clusters or the metallic- 
ity biases in red or blue giant stars), as well as uncertainties 
in determining uniquely their kinematics. Stars in the solar 
neighbo urhood have well-determined kinematics (e.g. GCS, 
RAVE: iNordstrom et al.ll2004l : ISteinmetzHJOOsh . however, 
they contain hardly any halo stars prohibiting an uncontam- 
inated and unbiased halo sample (Kinman 1995). To explore 
the halo in si t jj, intrinsically brighter objects are needed 
l|Sandagelll97Q) . but for those mainly line-of-sight velocities 
are available and many distance schemes associated with 
these stars present metallicity dependent systematics which 
can p roduce a false gradient in ki nematics across metallicity 
bins (|Fermani fc Schonrich 2013l'l. Among the most popular 



tracers are BHB stars (ISirko et al. 20041). RR Lv rae stars 



and globular clusters fe.g. lOdenkirchen et al.ll200H ). 

Distances for these stars are obtained spectro- 
photometrically: this implies that the uncertainty on the 
distance correlates with how noisy is the spectrum of a star. 
To assign an absolute mag nitude to the star one needs to 
identify its luminosity type IjBeers et aLlbOOd ) and therefore 
to have an accurate assessment of its surface gravity log p. 
This danger of misclassification is especially an issue for the 
noisy spectra of remote and hence faint objects deeper in 
the halo (e.g. lYannv et aLll2000l : ISirko et al.l 120041 '). For ex- 
ample, parameters in automated pipelines can be drawn into 
unphysical values or into points that are a local minimum 
by averaging over the mo st likely points, in case the pdf is 
not unimodal (see ^4.4.3. JLee et al.ll2008al '). 

It is therefore imperative to have a complete under- 
standing of how such uncertainties propagate in our anal- 
ysis and thus to bring possible biases under control: re- 
markable efforts have been made in this direction, but a 
complete set of rigoro us met h odolo gies has not been com- 
piled yet. For example iRvanI (|l992l ) demonstrated that er- 
rors in proper motion due to overestimates in distance can 
easily bias the kinematic measurements towards a c ounter- 
rotating trend. ISchonrich. Binnev fc Asplundl ([2013), here- 
after SBA12, lift the need to assume the shape for the 
velocity ellipsoid when estimating the distance of solar- 
neighbourhood stars by combining proper motion and line- 
of-sight velocity (iiioa): they achieve this by exploiting corre- 
lations between the measured U, V, W components of veloc- 
ity introduced by syste matic distance errors and iteratively 
correc ting the distance. ISchonrich. Asplund fc Casagrande I 
11201 If) showed that the c laim of a retrograde outer-halo in 
ICarollo et al.l (|2007l . [2010|) was a consequence of inappropri- 
ate treatment of distance errors. 

The aim of this paper is to achieve similar clarity in 
the context of BHB halo stars where only uios and tentative 
proper-motion information are available. Optimal control of 
systematics is essential to obtain a meaningful result in this 
sample. We attempt to meet this goal by contrasting two 
independent sample selections and four different estimators 
for halo rotation. 



In Section 2 we introduce this set of rotation estima- 
tors, and assess the reliability of each method depending on 
whether they use Vios information only or the full 3D mo- 
tion. Section 3 describes the two selections of BHB stars that 
we use as tracers of the stellar halo, and Section 4 presents 
the application of the rotation estimators on these samples. 
Section 5 compares our results with previous works, while 
Section 6 critiques the sample of Xll in detail. We sum up 
and look to the future in the Conclusions. 



2 METHODOLOGY 

We critique the model used by Dll. Assuming a radial power 
law potential and a distribution function Dll estimates the 
rotation and anisotropy of the stellar halo from a sample 
of BHB field halo stars for which only vios is available by 
marginalizing over proper motions. As an alternative, we 
assess the possibility of directly estimating rotation using ei- 
ther nios only but also full 3D motion. The latter estimators 
extract the average rotation and therefore are suitable even 
at large proper motion errors, as long as the systematic bias 
is small. We conclude the section by recalling general princi- 
ples of consistency for the rotational signature of a smooth 
component: these will later be useful to test the results of 
our analysis. 



2.1 Likelihood analysis 

One way to extract morphological and kinematic informa- 
tion from a sample is to fit all the data to a sufficiently 
general theoretical model. Commonly studied parameters of 
the Milky Way stellar halo are the rotation and the deviation 
of the velocity ellipsoid from a spherical shape, usually with 
the radial velocity dispersion varied against the two other di- 
rections. A generally accepted method is to assume a density 
profile for the stellar component and an overall potential for 
the Galaxy. Then one builds up a phase-space distribution 
function (DF) in the steady-state approximation, which will 
depend on parameters identifiable with anisotropy and ro- 
tation. One either can marginalise over observables that are 
considered too noisy or, more properly, convolve the distri- 
bution function with an error function replicating the obser- 
vation process. To estimate a value for the physically mean- 
ingful parameters, one can look at the likelihood of the data 
given the model, varying the free parameters. The results 
must be robust to the uncertainty in the data and give con- 
sistent results for different subsets, and the model needs to 
be flexible enough to fit various behaviours and trends. 

Dll also used BHB stars drawn from SDSS (DR7) to 
study the kinematics of the Milky Way stellar halo: thus, we 
adopt their model first to ensure consistency of our results. 



^ By in situ we refer to stars that are actually observed there, 
rather than stars which are extrapolated to be there via orbit- 
integration. 



2.1.1 Methodology of Deason et al. 2011 

Dll assume a radial potential of the type "I> ~ r~'' , normal- 
ized so that the escape velocity at the Sun is the one esti- 
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mated by ISmith et all (|2007l fl The density profile of BHB 
stars is assumed to be also of power-law type: p ~ r^°. 

To estimate the rotational signature and anisotropy 
of the sample, they adopt a double power-law distribution 
function depending on two parameters: /3 (an isotropy pa- 
rameter as defined in iBinnev fc T romainc 2008, §4.3.2) and 
7] (rotation), such that r] = corresponds to complete pro- 
rotation and r; = 2 to complete counter-rotation: 

F„,^{E,L,L.) = c^ (l + (1 - r?)tanh (:^)) L'^^ E^ , (1) 

where E is the binding energy, A is a smoothing parame- 
ter, C/3 is a constant that enables us to make the real-space 
density independent of (3 and 



C^ 



Pil - 2) 



+ 



(2) 



Assuming dominant errors on the proper motions, the 
DF is marginalized over transverse velocities. The free pa- 
rameters (77, /3) are then estimated via Markov Chain Monte 
Carlo (MCMC) sampling on the parameter space with the 
probability density being the likelihood of the dat4j given 
the model, so the probability density in parameter space is 
the exponential of: 



L{Ti,P) = J2^ogF^,fi{k,h,d, 



ll f^lOS,l/5 



(3) 



where 



F^^I3{1, b, d, vios) = / dvidvtFr,^i3{l,b,d,vi,vi,,vios)- (4) 



See Appendix A for details on the convergence of the 
MCMC estimate and for a test of the algorithm implemen- 
tation. 



2.1.2 Notes on the method 

The choice of an index a ^ 3 for the density profile makes 
it impossible to compute a normalization constant for the 
DF in lU as it contains an infinite mass at the centre: 
J'd^a;d^ii Fr,^i3{E,L, Lz) = co. Therefore, in order to be 
able to compare different models we must fix the density in 
the sampled region of real space. Thus 



p(r)= / d\F^AE,L,L,) 



(5) 



is independent of /3. 

A convenient consequence of fixing the density is that 
we can neglect selection effects which would otherwise apply 



^ The Sun is taken to be at a distance of 8. 5 kpc from th e Galac tic 
Center. The escape velocity reported by ISmith et al.l 1120071 ) is 
''osc = 498 — 608 km s^^ with 90% confidence and an average of 
544 km s~^ . We adopt the average value. 

^ Heliocentric galactic polar coordinates are converted into galac- 
tocentric cylindrical coordinates reference frame, by assuming 
(x,y,z)Q = (8.5, 0, 0) kpc, Vcq = 220 km s~^. We are aware 
that this specific parameter combination is at odds with the mea- 
sured proper motion of Sgr A* , but keep those values consistent 
with previous studies. We further assume a solar Local Standard 
of Rest velocity of ([/, V", iy)(^=(ll. 1,12 24,7. 25) km s'^, as up- 
dated bv lSchonrich. Binnev fc DehnenI 1 I2OIOI I. 



JMcMillan fc BinnevI [20111 ). All the same, we are left with 
the uncertainty on the real density profile of the tracers. If 
the steepness of the profile is overestimated, for example, 
the model considers distant stars less probable than they 
really are. This results in these objects having more weight 
than their closer companions. An unfortunate aggravation 
of this bias is that fainter/more remote stars are known to 
be noisier and are intrinsically a bad tracer of the rotational 
signature. 

When marginalizing over proper motions with /3 > 0, 
we compute an integral that involves a singularity in the DF 
IT]) at L = and the treatment of this singularity is crucial 
to retrieve /3 correctly. 

Fortunately, the singularity is integrable for /3 < i. 

The estimate of the anisotropy is still a matter of con- 
cern though: /3 appears as L~'^^ . This term gives rise to 
a bimodality in velocity space, which cannot adequately 
fit the single-peaked velocity distribution observed for lo- 
cal halo members. Hence any anisotropy measurement with 
this model is condemned to fail, in most cases with an az- 
imuthal bias caused by depopulating the velocity distribu- 
tion around V ^ 0, which also explains why the anisotropy 
is at odds with observations of a strongly radially dominated 
velocity e llipsoid for local stars (see e.g. ISmith et al.l 120071 : 
[Schonrich. Asplund fc Casagrande Il201ll 'l. 

Equation (6) of Dll states that the mean rotational 
velocity for model IT} depends on radius only. By computing 



(v^) 



Jd'^vv^Fr,,p{E,L,L,) 
fd\Fr,AE,L,L,) 



(6) 



we see that this is true only in the limit A — >■ 0. Indeed, 
one can show that if tanh(I/z/A) 7^ sgn(Lz), then (v^) de- 
pends on the full 3D position. We will therefore adopt this 
assumption even though it was not stated in Dll. 

Furthermore, we note that the choice of the potential- 
density pair gives rise to an unphysical mean velocity in the 
inner Galaxy. According to eq. 6: 



hm \{v4,)\ = cx). 

7'— )-0 



2.2 Model-independent estimators of rotation 

We present here three model-independent estimators of rota- 
tion: one relies on wios informati on only (a similar app roach 
has been previously adopted bv lFrenk fc Whitelll980l ) and 
the other two on the full 3D motion. SDSS proper-motion 
at these distances are of the size of the noise associated with 
them. All the same, given that the estimators rely on aver- 
ages on the sky, valuable information can still be extracted 
from them provided that their errors are random and not 
systematic. To set ourselves in such a condition we correct 
for the known s ystematics of astrometric frame-dragging 
(jSchonrichl 120121 ); if there is any other significant system- 
atics left, the estimators will be biased differently and will 
therefore contradict each other (see SBA12). Further, one 
of them outputs the U-component of Sun velocity and the 
latter information can be directly compared to the values in 
literature to assess the presence of proper-motion systemat- 
ics significant enough to bias the estimator. 
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2.2.1 iiios estimator 

Inferring the rotational signature using only wioa presupposes 
knowledge of the global velocity distribution: i.e. how the 
velocity ellipsoid projects into wios. In fact, it is sufficient to 
assume that the azimuthal velocity is not correlated with 
the other components (the velocity ellipsoid tilts between U 
and W , but not V) and that the means of the other two 
components are ~ 0. This minimal set of assumption allows 
us to measure rotation from the de-projected vios on the 
azimuthal direction: 



^loa = (V - V0) • S => V • S = Wlos + V0 • S, 



(7) 



where s is the l.o.s. unit vector, e^ is the unit vector in the 
azimuthal direction of Galactic cylindrical coordinates and 
V0 is the Sun's velocity w.r.t. the LSR. Therefore: 



iv^) 



«los + Vq ■ S 



(8) 



Hereafter, we shall refer to the de-projected l.o.s. veloc- 
ity as 



71- («los) 



Wloa -I- Vq ■ § 



If v^ is independent of position we can take a global av- 
erage and expect that V^{v^) ~ 0, where x is the 3D posi- 
tion vector. The best estimate of the mean streaming motion 
is then obtained by averaging the deprojected measurements 
scaled by their associated variances following a maximum 
likelihood argument. In practice, we fit a straight line in the 
Xj, ■K~^{v\os) plane with the "weighted" least square routine 
of gnuplot having marginalized over the other two position 
variables. If 'K~^{vias) depends on position, the estimator 
still gives a best estimate of the position-dependent rota- 
tion velocity. The variances must account for the projection 
of the U , W components of velocity along the l.o.s. as well 
as for the dispersion in the V distribution av, so that the 
final variance is: 



^C 






(9) 



where cr^j^^.j is the radial velocity measurement error for the 
C,-th star, and the sum describes the projection of the veloc- 
ity ellipsoid into the l.o.s. with direction vectors e^ and dis- 
persion Oj . The deprojection factor has a singularity, where s 
is perpendicular to e^. In Fig[l] we show s ■ e^ in the (x, y)- 
galactocentric Cartesian plane at altitudes of \z\ — 5 and 
\z\ = 20kpc. A deprojection factor of 0.1 means that a rota- 
tion velocity of 20kms~^ leaves a signal of just 2kms~^ in 
the mean l.o.s. velocity. Hence the region where l.o.s. veloc- 
ities can be reliably used is confined to small galactocentric 
radii and low altitudes with the best regions along the y-axis. 
The variation also implies that weighting the data points 
according to their uncertainty is mandatory to obtain any 
meaningful results. Neglect of a^ will lead to large errors (see 
i|2.3.1|) and fiuctuations of the results|j Fortunately we do 
not require exact knowledge of the uncertainty. Since varia- 
tions in the numerator of (|9} are comparably small and since 





depr. fact. 




* We suspect this is the case in study of iHattori et al.l 1(2013), 
who do not mention any care towards this issue. 



depr. fact. 



Figure 1. Colour coded deprojection factor in the {x,y)- 
galactocentric Cartesian plane at altitudes of \z\ = 5kpc (top 
panel) and \z\ = 20kpc (bottom panel), limiting W x^ + y'^ < 
45kpc. The white cross denotes the Sun's position; note the dif- 
ferent colour scales of the plots. 



weighted least squares only improve the statistics but do not 
introduce any bias, we can safely set ^■^^{o'jej ■ s)^ con- 
stant, as long as we recover the approximate shape around 
the singularity. 



2.2.2 3D estimators 

We have extracted rotational information out of wios alone 
both via a model and via a direct estimator. i;ios does not 
hold the full kinematics though and we could be missing 
stream-like structures in our sample with peculiar motions. 
Therefore to disentangle the possible bias by streams, we 
need to make use of the full 3D velocity. 

This implies that the major part of velocity information 
actually comes from proper motions. While velocities of in- 
dividual stars are not reliable - the error is approximately 
as large as the signal - large number statistics solve that 
problem for aggregate quantities like the average rotation|j 
The error in the mean quantities is anyway best estimated 
from the as-observed scatter and does not require a-priori 
knowledge of the proper motion (random) errors. 

More of a concern are systematic errors on the proper 
motions like astrometric "frame-dragging" . We checked that 
neither the corrections from iSchonrich (,2012 ) nor differences 



^ Even if the noise is actually bigger than the signal, which is 
the case for SDSS measurements of stars beyond 20 kpc from the 
Sun, this introduces no preferred direction of motion and hence 
no bias. 
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between the data releases can alter our results, for details 
please see subsection 14. 1.1 1 

Furthermore, SBA12 show that the mean streaming mo- 
tion can be recovered from the position-dependent compo- 
nent of velocity (U) in the local frame: 



U 



vr cos a + v^ sin a, 



(10) 



where the velocity vector of the star v is projected by e^, 
onto the local radial direction, u^ is the effective rotation 
and 

/ d sin Z cos 6 \ ,,,, 

a — arctan ; ; (11) 

\Ro — a cos ( cos b J 

is the angle between the projection onto the plane of the 
long axis of the velocity ellipsoid and the Sun-GC line (see 
equations 26-28 in SBA12). The average value of U depends 
on sin a, so that we can estimate the global mean w^ from 
the linear regression: 



IJJ) = (?;0)sinQ - [/q. 



(12) 



The linear fit recovers both the mean rotation and the 
[/-component of Sun's velocity, [/q- 

Finally the direct estimate of the mean streaming ve- 
locity of the sample can serve as consistency check that the 
[/-estimator is not biased from proper-motion systematics: 
if it were, contradictory results would emerge between the 
two estimators given their different bias. 

2.3 Sample systematics 

A broken kinematic signature across different sub-divisions 
of a sample is the consequence of either a distance error, 
a pipeline error or a stream-like structure. Provided that 
tracers of all prominent unrelaxed sub-structures have been 
removed from the sample, no matter how we sub-divide our 
sample, each sub-set must have the same kinematics. What- 
ever portion of the sky we look at, we have to obtain the 
same rotation estimate, within the uncertainty. In partic- 
ular, splitting the sample in random sub-samples decreases 
the statistical signal-to-noise and consequently enhances the 
allowance for discrepancy in the signature between two bins. 
Therefore, an anomaly among different sub-divisions of the 
sample is even more concerning and serves as a strong indi- 
cation for a bias in the analysis. 

2,. 3.1 Distance errors 

When we recover the rotational signature of a population 
of stars where only wios is available or reliable, an error in 
distance induces changes in both the angle between the line 
of sight and the vector from the star to the Galactic Cen- 
tre, and the projection of wios onto the plane tangential to 
the direction to the Galactic centre. Thus a distance error 
influences directly the rotational signature. This will not be 
a constant systematic, but a function of the position on the 
sky. Distance errors shift stars en masse in velocity space. 
Distance over-estimates tend to shift velocities upwards to- 
wards the escape velocity, thus narrowing the range of phys- 
ically accessible velocities to be considered. An increase in 
velocity will nearly always reduce the likelihood of a star 
when a model is being used to calculate the probability of 
the data; by carrying a star beyond the escape velocity they 
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Figure 2. We move each star in the FS12 sample along the l.o.s. 
to put it at distance d, then calculate {l/(s • e^)} and plot it ver- 
sus d (green line). The red line represents the minimum l/(s-e^): 
both are almost perfectly linear functions of distance beyond 10 
kpc. Blue crosses depict the stars from FS12 at their actual he- 
liocentric distance. For graphical convenience we do not show the 
6% of stars which exceed the vertical range. For better compari- 
son (see §5) with the analysis of Deason et al. (2011), we extend 
the distance range to 90 kpc. 



can even reduce the star's probability to zero. When a direct 
estimator is used, a distance over-estimate tends to increase 
the computational weight of a given star. 

This effect does not change the sign of the rotation con- 
tribution of individual stars, but it distorts the values with 
a significance that scales with distance from the observer. 
In Fig. [2] we see that the factor by which the star's veloc- 
ity is multiplied increases almost monotonically heliocen- 
tric distance. Indeed, distant stars are intrinsically a bad 
tracer if proper-motion is not available: at large radii the 
streaming motion v^ contributes very little to «ic,s; in fact 
wioa -|- V0 ■ s ~ n,, and v^ is dominated by the tangential 
component of velocity (hereafter v^). In order to change 
wios by a small amount, v^ needs to change significantly: the 
signal-to-noise scales inversely with radius making rotation 
estimated via vi^s of distant stars, dominated by noise. 



2.3.2 Stream-like structures 

A discontinuous rotational signature can hint at the presence 
of a sub-structure in the data: a stream for example can lead 
to inconsistencies between different rotation estimators, esp. 
for the vioa based estimator and model. Even though we do 
not hunt for substructures here, we must allow for uniden- 
tified accretion remnants mimicking rotation. To break the 
degeneracy between biased data and a real stream, we need 
to identify what is characteristic of a substructure and what 
is not. An unrelaxed substructure is identified as such by 
its peculiar signature in parameter space (both physical and 
kinematic) . Correlations between kinematics and metallicity 
or between colour and kinematics are examples of signatures 
of a possible accretion event. By contrast binning in lumi- 
nosity type acts like a random selection w.r.t. kinematics 
and hence a coherent signature across the bins must be ob- 
served. Further, provided that all detectable over-densities 
in observable space have been removed, binning the sample 
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spatially implies that even if there are remnants of accretion 
events left they are washed out by the smooth component 
and thus a coherent signature is expected across all the sub- 
samples. 



2.3.3 Pipeline Systematics 

A bias in the template fitting of the pipeline can break the 
rotational signature as we move from the left (/ < 180) to the 
right {I > 180) side of the Galaxy. The spatial dependence 
will make a genuine coherent streaming motion stronger on 
one side and faint, or even of opposite sign, on the other 
side. In general, a higher degree of caution must be adopted 
especially for hot metal-poor stars, since their spectra are 
almost free of metal-lines and a line-of-sight velocity deter- 
mination based on the very broad Balmer lines is vulnerable 
to template misfits. 



selection process (see i|3.2[) . To both samples we apply geo- 
metric cuts to exclude the disk and Sagittarius. 

Astrometric accuracy is of crucial interest in our study: 
the statistical error on SPSS proper- mot ions is of the or - 
der of ~ Smasyr"^ (|Munn et al.ll2004l H: [Pong et al.ll201ll '). 
We constructed our statistics (see i|2.2.2l) such that we do 
not require any knowledge on the proper-motion error distri- 
bution: the only implied assumption in our error estimates 
is that the observed velocity distributions including the er- 
rors follow a roughly Gaussian distribution. The violation 
of this assumption does not affect our parameter estimates, 
but leads to a mild misjudgement of their errors. Proper- 
motion systematics are more of a concern: we assess their 
effect in our analysis in detail in i]4.1.1l 



3 SAMPLE SELECTION AND CALIBRATION 

To study the Galactic halo in situ we need tracers that are 
bright enough to allow for distances > 5 kpc. The studied 
stars should have relatively well-determined absolute mag- 
nitudes and should not be easily mixed up with other types 
of stars, e.g. foreground dwarfs from the Galactic disc(s). 
Late type giants are problematic by their easy confusion 
with dwarf stars and by their steep relation between colours 
and magnitude, which also amplifies reddening uncertain- 
ties. The opti mal objects would be RR Lyrae stars (cf. 
iKlemenj |2010| V but in addition to spectroscopy they re- 
quire a good coverage in time to resolve their oscillations 
and hence are not globally available (Sesar et al.ii20 10V 

Blue horizontal branch (BHB) stars have the advan- 
tage of a relatively narrow absolute visible magnitude dis- 
tribution that shows only weak dependence on metallicity. 
It should be kept in mind, however, that there is the blue 
hook with fainter objects on the hot side of the BHB, and 
that, despite most Galactic stars being redder in colour, they 
mingle in observations with the fainter blue stragglers and 
main-sequence objects. Three techniques have been devel- 
oped to identify B HB stars by looking at hydrogen lines: 
meas uring the size l|Philip fc Adelmanlll994l ')) or the steep- 
ness (|Chalongei Il973l ) of th e Balmer Jump and us ing the 
width of the Balmer lines l|Searle fc Rodgerg|l966l ). After 
identification BHB stars are treated as standard candles. It 
should be kept in mind, however, that the use of BHB stars 
biases any sample to old, metal-poor populations. Moreover, 
theoretical models show a metallicity dependence of order 
~ 0.3 mag in the V-band magnitude (and similarly in the 
Johnson g-magnitude) , emphasizing the need to control any 
result not only for co ntamination, but also for abs olute mag- 
nitude uncertainties IjFermani fc SchonrichlbOlSl '). 

We construct two samples from SDSS DR9: one selected 
via filters in colour, surface gravity and effective temperature 
and another selected by Xll via direct use of Balmer lines. 
We will refer to the first sample as photometrically selected 
calling it FS12 and to the second sample as spectroscopi- 
cally selected, calling it Xll. Even though the FS12 sample 
is drawn from data for which spectra are in fact available, 
the chosen labels aim at emphasizing the difference between 
using and not-using Balmer-line shape measurements in the 



3.1 Selection via photometry 

We c onstruct a samp le of BHB stars, drawn from SDSS 
fDR9 lAhn et al.ll201^ i p. The selecti o n ben efits f rom the ex- 
tensiv e work done by lYannv et al.l (|2000|) and ISirko et al.l 
(|200J), who revised the BHB identification procedure work- 
ing with SDSS data. 

Both agree on an initial colour filter such as: g — r £ 
[—0.4,0] and u ~ g £ [0.8, 1.4]. The comparison with ultra- 
violet colour u — g provides indirect information on the 
Balmer jump and hence is an indirect gravity estimator. Still 
in this colour regime (which we will denote TV) there are con- 
taminations from other A-type stars (mostly blue stragglers) 
and from F stars due to "intr insic variations and photomet- 
ric errors" (jSirko et al.l 120041 ). While F stars can be filtered 
out by their lower temperature, A-type stars mostly differ by 
their surface gravity, so that both groups consequently use 
the Balmer line width to separate them. Their understand- 
ing of the problem can be summarized as follows: main- 
sequence A stars can be distinguished from lower surface- 
gravity BHB stars (which lie in the region log ff € [2.5, 3.2]) 
by lo oking at the width of the HS lines alon e ( Yannv et al.l 
liooO) or H5 and H-y lines (|Sirko et al.ll2004l '). The discrimi- 
nation is not sharp and one must still expect contamination 
of order ~ 10% (jSirko et al. 2004). The situation deterio- 
rates towar d s fain ter magnitudes and hence noisier spectra: 
ISirko et al.l (|2004l) estimate contamination of order ~ 25% 
at p > 18. Indeed iLee et al.l (|2008b|) show that the low spec- 
tral resolution in SDSS and limited signal to noise lead to 
larger uncertainties in stellar parameters including log{g) 
up to a level where subgiants and dwarfs get mixed up. This 
represents a serious challenge for the proposition that stars 
can be cleanly sorted in type using only the log(.(?) estimates 
from the SEGUE s tellar pipeline (e.g. ICarollo et aLllJOOTJ , 
I2OI0I : iDeason et allfeoil ') . 

According to the dangers reviewed in the previous para- 
graph, we minimise the risk of contamination from other 
luminosity types, by selecting: 



^ See Appendix B for the SQL query. 
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g<i8 

2<log5< 3.3 

0.8 < {u~g)o < 1.4 
-0.4 < (g - r)o < 
7250 < Teff/K < 9700 
\z\ > 4kpc & r > lOkpc 
(a,5)^Sgr 



(13) 



where the last two criteria are geometric cuts to exclude disk 
stars and members of the Sagittarius stream (e.g. Dlljj. 
The estimates of surface gravity, ef fective temperature and 
metal licity are from the analysis of IWilhelm. Beers fc Gravl 
l| 19991 ) as it was specifically designed for hot stars. This leads 
to a sample of 1585 stars. In Fig. [3] we plot the geome- 
try of the sample: the top panel shows the distribution in 
the {x, 2/)-plane colour-coded according to the z-coordinate 
(Cartesian Galactocentric reference frame) and in the bot- 
tom panels we show the heliocentric (left) and galactocentric 
(right) distance distributions. 

We further investigate the purity of our sample in the 
colour-colour plane and find grounds to suspect the presence 
of the blue hook in the stripe —0.3 < (g — r)o < —0.2 (see 
Fig. 21 which suggests that a,u — g dependent cut in the g~r 
colour might more cleanly identify the horizontal branch. To 
a zero-th order approximation, a constant cut at {u — g)o > 
1.15 (rather than {u — g)o > 0.8) is cost-effective and will 
therefore be our initial choice to calibrate the reliability of 
H13p . The tighter cut in {u — g)o shrinks the sample to 1297 
objects. 

Even though our stringent selection criteria cause a 
drop in the sample size w.r.t. previous works (e.g. the sam- 
ple of Dll has ~ 3500 stars due to more generous filters in 
surface gravity and apparent magnitude) , none of our filters 
correlate with kinematics and thus they do not prejudice our 
analysis. They also do not bias the distributions of physical 
parameters such as for example metallicity: the distribution 
peaks at [Fe/H] = —1.9, in good agre emen t with the esti- 
mates from the samples of IXue et al.l (|2008l . l201ll ). who do 
not impose any cut in surface gravity or apparent magnitude 
for example. Further, we judge the increase in the statisti- 
cal noise in the model fits a convenient price to pay in light 
of the reduced risk of contamination. The fact that our se- 
lection criteria in logg and p-band magnitude do not bias 
the inferred kinematics is confirmed by a second sample se- 
lected via more generous photometric filters, but combined 
with Balmer-line shape measurements, which presents an 
homologous rotational signature (see following Section). 
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Figure 3. Geometry of the FS12 sample (top panel): distribution 
in the (x, j/)-plane with the origin at the Galactic Centre, altitude 
z is colour-coded (Cartesian Galactocentric reference frame). 207 
out of 1585 stars (13%) are in the southern galactic hemisphere. 
The bottom panels show the heliocentric (left) and Galactocentric 
(right) distance distribution. The dark blue histograms refer to 
the FS12 selection I I13I I. while the light blue ones show the sample 
without the geometric cuts. 
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Figure 4. BHB stars drawn from SDSS DR9 included by the 
selection I I13I1 in the («— 9)0, {g—r)o plane. Every point is coloured 
according to the metallicity of each star. The commonly adopted 
colour cuts appear to include some part of the blue hook. 



3.2 Spectroscopic selection 

XI 1 presented a spectroscopically selected sample of ~ 5000 
BHB stars drawn from SDSS DR8 and which has been 
kindly provided to us. The advantage w.r.t. a selection via 
photometric information is that the three independent meth- 
ods of identification separate BHB stars from blue strag- 
glers and main-sequence stars more cleanly by combining 
colour cuts with Balmer-line shape measurements: the scale 
width-shape method fits Sersic profiles to the Hy lines (cf. 

'^ The masking of Sgr is performed in right ascension-declination 
space (a, 5) according to the polygon kindly provided by the au- 
thors of Dll on request (see Appendix B). 



IClewlev et a"Lll2002l : ISirko et al.ll2004 IXue et al.ll2008l . I2OIII ) . 
The parameters that describe such a profile are essentially 
the amplitude (a-y), the dispersion (6^.) and the steepness 
(c-y) of the modified Gaussian that is fitted to each H^ line: 



y{\) = 1.0 — a-y exp 



lA-Aol 



(14) 



We retrieve Xll's stars in DR9 and use the astrometry 
of the latest data release; we exclude contributions from the 
disk and the region of Sagittarius by applying the same geo- 
metric criteria of the previous section: this leads to a sample 
of 2563 stars (2455 with [Fe/H] £ [-3, 0]). We show its geom- 
etry in Fig. [5] Among these ~ 1400 have c-y ^ 0.95, which 



© 2013 RAS, MNRAS OOO.[TH20l 



8 Fermani & Schonrich 




20^ , ,40 60 
d / kDC 



20 ,,40 60 
r / kDC 



Figure 5. Same as FigO but for the Xll sample. Here 526 out 
of 2563 stars (21%) belong to the souther galactic hemisphere. 



is the limit value for which the method of [Clewlev et al 



2002 ) has been calibrated. Given that neither Sirko e t al. 



2004 ) nor lXue et al.l l|201ll ') have re-calibrated the method 
in the region c^ ^ 0.95 we cautiously flag these stars and 
first consider them in our analysis but later analyse them 
separately. 

We retrieve 1400 stars of our photometrically selected 
sample in the sample of Xll, which implies that almost 90% 
of the FS12 sample is included in Xll: the latter sample 
was selected using a classification based on both Balmer- 
lines and colour cuts and therefore constitutes the more ro- 
bust identification of BHB stars to the current understand- 
ing. Cons equently, c r oss-m atching our FS12 sample with the 
sample of lXue et al.l (|2011h provides a rough estimate for the 
level of contamination from non-BHB stars that we can ex- 
pect: i.e. ~ 1 0%, which very we ll matches the independent 
prediction of ISirko et al.l (1200J). Note that when we com- 
pare the sample of Dll with Xll, only 47% of their stars 
have Balmer-line shapes identifiable with BHB type stars, 
even though Dll drew their sample from DR7 and Xll is 
based on DR8, which is significantly larger. If we exclude 
high surface gravity and faint stars from t he sample of Dll 
in DR7, the proportion of stars matched bv lXue et al.l (|201ll ) 
increases to 71%, meaning that the maj ority of the stars we 
rejected are non-BHB stars according to lXue et al.l pOlll lrl 



* Note that the value reported in Dll at the beginning of their 
§3.1 in these regards can be s omehow mislead ing. The authors 
report that 88% of the stars of |Xue et all |200S| ) (based on DR6) 
is found in their larger sample. But this actually mean s that only 
41% o f the stars in Dll sample is found in the one of IXue et al.l 
(I2OO3) and we do not have any reassurance on the level of con- 
tamination for the remaining 59% of the stars in their sample of 
(~ 3274 objects). 



3.3 Absolute Magnitude 

The second issue that needs careful handling is the es- 
timate of the absolute magnitude for the BHB stars. 
ISirko et al.l (J200J) show that the assumption of a con- 
stant luminosity and m ass, temperature-independent for 
Tcft G [8000, 12000]K, (JBaev. Markov fc Spasswal I2001D 
would lead to an average overestimate in distance modulus of 
0.18 mag with respect t o a temperature-dependent relation - 
ship for the luminosity (|Dorman. Rood fc 0'Connelllll993l ). 
Fu rthermore, the popul ar value of Mg = 0.7 was calculated 
by iLavden et al.l l|l996l ) for halo RR Lyrae stars and then 
adopted for BHB stars given the similarity in their absolute 
magnitude. The justification of this app roximation depend s 
on the blue re s ponse of the filter in use (jYannv et al.ll200Gl ). 
IPreston et al.l ([199l|) report that BHB stars in the blue end 
of TZ are 0.7 mag fainter in g than the RR Lyrae star^. 
Therefore the uncertainty reported for the absolute magni- 
tude of BHB stars needs to be interpreted generously given 
the non-standardization of the blue filters and the conse- 
quent difficulty in comparing estimates from different sur- 
veys. 



iFermani fc SchonrichI l|2013l ) show that the assumption 
of a constant Mg for BHB stars introduces systematics cor- 
related with colour and metallicity and propose a colour and 
metallicity dependent distance calibration for BHB field halo 
stars: 

Mg{ [g - r)o, [Fe/H]) = 0.0075 exp(-14.0(<7 - r)o) + 

-H0.04([Fe/H]+ 3.5)^ + 0.25. (15) 

This formula covers the theoretical expect a tions from the 
BASTI isochrones l|Pietrinferni et all [20oi , |2006| ). and in 
statistical tests on field stars proves t o be more accu- 
rate than other avail able calibrations (e.g. ISirko et al.ll2004l : 
iDeason et al.M2011bl ). which were only based on globular 
clusters. In particular, (|15|l hits the precision threshold al- 
lowed by the sample size, on FS12 the average fractional 
distance error is 0.01 ± 0.04 and 0.02 ± 0.03 on Xll, and is 
the only calibration to pass the "falling sky" test of SBA12 
(jFermani fc Schonrichll2013l . §4.1 and 5.6 respectively). Cru- 
cially in this context, it covers the metallicity dependence 
and hence avoids overestimating the distance to metal-rich 
objects and underestimating that to metal-poor ones result- 
ing in metallicity-dependent kinematic biases. We will adopt 
the above formula throughout the paper, unless otherwise 
stated or when comparing to Dll, who assume a constant 
BHB absolute magnitude. 



4 ROTATIONAL SIGNATURE OF THE MILKY 
WAY STELLAR HALO 

We present below four independent estimates of the rotation 
of the Milky Way stellar halo, traced with the two different 
samples of BHB stars presented in the previous section. 
The likelihood analysis yields a very weak rotation on 



^ Both lLavden et al.l 1119961) and lPreston et al.l |l99ll') report re- 
sults for the V magnitude, which lYannv et al. 1120001 ) usefully 
converts into the q m agnitude via the colour transformation of 
iFukugita et al.l l ll996l) . 
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all 


[Fc/H] e [-3, -2] 


[Fo/H] e [-2, 0] 


FS12 

(^-M^los)) 
{Vtl,)u-eBt 


0.83 ±0.08 
17±9kms-i 
0±llkms-i 
6±6kms-i 


0.84 ±0.17 
ll±15kms-i 
2 ±17 kms" 1 
-3±10kms-l 


0.79 ±0.12 

22±12kms-l 

-l±14kms-l 

13±8kms-^ 


Xll 

(^-^(^los)) 
{v^)u-eBt 


1.04 ±0.08 

-3±7kms-i 

-13±10kms-l 

-7±5kms-l 


1.22 ±0.15 

-18±llkms-i 

-7±16kms-i 

-4±8kms-i 


0.90 ±0.10 

9±9kms-l 

-18±13kms-i 

-9±7kms-i 



Table 1. Estimates of rotation using both i^ios information only and full 3D motion on two samples drawn from SDSS DR,9. r] is the 
rotation parameter from model JlJ and implies prograde rotation for »? < 1 and retrograde motion for rj > 1. {'t~^{vi^b)) is the rotation 
estimate from the i)iog-estimator and (i'rf,)[/_est the one from the [/-estimator (see §2.2). 



both samples, consistent with non-rotation. On FS12 we do 
not observe any significant trend in the rotational parameter 
rj (see ^2.1.1|l across metallicity: ?7[Fe/H]<-2 - '7[Fo/h]>-2 = 
0.05 ± 0.21, while there is a gradient in the Xll sample: 
'?[Fe/H]<-2 — '7[Fo/H]>-2 ~ 0.32 ± 0.18. This mismatch dis- 
appears on the sub-sample where FS12 and Xll overlap 
(~ 1400 stars). The uios estimator confirms all the above 
conclusions. Fig. [6] presents these results together with the 
anisotropy estimates from the model in the different cases: 
FS12 favours an anisotropy consistent with zero, while Xll 
suggests a very weak radial bias: /3 — 0.08 ±0.05. It is worth 
remarking that as pointed out in §2.1.2, the model is in- 
trinsically inadequate to estimate the anisotropy due to the 
L~^^ term - therefore, the obvious discrepancy to the ra- 
dial bias in local samples with full 3D information is just a 
sign for its failure. We note that adopting either a constant 
distance scale {Mg — 0.7, to be consistent with Dll) or the 
colour- metallicity dependent calibration (|15p does not alter 
the above conclusions. 

When we apply the U and v^-SD estimators, they con- 
verge to the same conclusion of a rotational signal consis- 
tent with zero. On one side, the fact that two 3D estimators 
agree implies that whatever the nature of the errors affect- 
ing the proper-motion, they are either random or not sig- 
nificant enough to bias the average the estimators rely on. 
Indeed, were there still significant systematics, the estima- 
tors would diverge. Further, the [/-estimator provides also 
the t/-component of the Solar velocity. On FS12 we measure 
Uo = 15.7 ± 7.4 kms"^ and Uo = 7.8 ± 6.8kms~^ on Xll: 
within the error bars both agree with the accepted value of 
Uo, again supporting the reliability of the 3D estimators. All 
estimat es are obtained applying the proper-motion correc- 
tion bv lSchonrichI (|2012l ): in i|4.1.1l we quantify the effect of 
neglecting this amendment. 

On FS12 all estimators agree in measuring no or very 
weak rotation in both metallicity subsamples (see Table [T|. 
However, on the Xll sample the ID-estimators measure a 
slight, though significant, discrepancy between the kinemat- 
ics of metal-poor and metal-rich stara^^l. 

For the entire halo population, all estimators measure 



^^ We call metal-poor stars the ones with [Fo/H] g [—3, —2] and 
metal-rich the ones with [Fe/H] G [-2,0]. The ratio is 9:20 for 
both samples. 



a rotational signature consistent with zero on both samples. 
On FS12, both the ID and the 3D estimators agree that 
there is no significant discrepancy in the rotation of metal- 
poor objects w.r.t. metal-rich ones. On Xll the situation is 
less clear as the ID estimators support the existence of a 
gradient in rotation across metallicity bins, while the 3D es- 
timators disprove it|^ This issue arises on the sample where 
we had already flagged 50% of the stars as potentially prob- 
lematic: in §6 we critique this discrepancy in detail. 



4.1 Notes on uncertainties 

4-1.1 Uncertainties from datasets 

Proper motions of stars have significant systematic errors, 
which may affect the 3D-estimators. The main systematic 
that affect the prope r motions are "frame-dr agging" and 
chromatic aberration IjKaczmarczik et al.ll2009l ). The astro- 
metric frame in Sloan is fixed using galaxies, and can be 
"dragged" by the inc lusion of stars i n the "galaxy" sam- 
ple (see dis cussion in ISchonrichll2012l ). The correction de- 
veloped by ISchonrichI ( 20121 ) is designed for usual stellar 



colours and the chromatic aberrati on for the blue r BHB 
stars may be different. However, the lSchneider et al.l (|2010r ) 
quasar sample harbours only very few objects with observed 
g ~ r < (about 10%) providing us with a meagre basis 
for tests. While for gi — r > the sample displays almost 
perfect stability, the blue end deviates at the 1 — 2cr level 
towards smaller corrections. Fortunately, the correction on 
the DR9 sample shifts the average rotational velocity (w^) 
by a moderate ~ 0.6— 1.6kms~^ for the (D^)-estimator and 
by ~ 1.1 — 1.6 kms~^ for the [/-estimator towards retrograde 
motion. Anyway, in most cases we are only interested in the 
rotation difference between two sub-samples, which is nearly 



^^ Note though that an aberration in the l.o.s. velocity would 
shift the estimate of the ID estimator leaving nearly unaffected 
the ones of the 3D estimators (see ^.l.ll l. An aberra t ion o f 
±5(10)kms~^ (following what measured by ISchonrichI 1120121 ) 
for blue metal poor stars, though for a different portion of the 
sky) will shift the •7r~'^(i)ios)-estimate of —18 ± llkms"^ by 
±9.7 (19.4) kms~^ while moving the 3D-estimator by only 1 (2) 
kms~^. 
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Figure 6. Rotation and anisotropy of BHB stars from SDSS 
DR9 estimated via likelihood of the model jTJ given the data. 
Green line is the one sigma contour for the full sample, red line 
is associated with the metal rich stars and the blue one with the 
most metal-poor objects. Full lines are associated to the sample 
selected via l|13|l . dashed lines refer to the spectroscopically se- 
lected sample from Xll. The model-independent t^iog estimator 
(see §2.2) agrees with the findings of the model. Both methods 
use nios information only. 



unaffected. Therefore, even with the proper-motion uncer- 
tainty the mean v^ provides valuable information to check 
the robustness of the model-dependent rotation estimates. 

The heuristic exercise of perturbing our wios by a sys- 
tematic ±5kms~^ shows that the SD-estimators are rather 
robust against systematics in the l.o.s. velocity: 



{'»<l>)vu 



(V^)v,, 



±lkms 



(16) 



The ID «ios-estimator is on the contrary affected by a shift 
larger than the actual systematic introduced: 



(tT ^(uios))„. 



,zb5 kms 



(«1os))d,. 



±9 kms" 



(17) 



This strong deviation might be surprising at the first glance, 
since the left and the right side of the sky should balance 
each other. However, the sample is strongly asymmetric, and 
worse, the geometric factors in the ID-estimator deproject- 
ing the observed mean motion into a rotational signal am- 
plify any systematic on the line-of-sight velocities. This con- 
firms the necessity of validating rotational estimates via SD- 
estimators and supports the conclusions of the latter when- 
ever they disagree with the ID-estimatorsr^l 

As a further check we repeated all experiments on 
earlier data releases DR7 and DR8, before SDSS un- 
der went significant changes espec i ally to the astrome- 
try (|Adelman-McCarthv et all l2008l : lAbazaian et all l2009l : 



^^ For example the kinematic trend across difi^erent metallicity 
bins found in the Xll sample and supported by the ID-estimators 
only will in fact result in being most probably an artefact (see §6). 





u- g > 0.8 
V 


/3 


u- g> 1.15 


n 


all 
rich 
poor 


0.83 ±0.08 
0.75 ±0.12 
0.97 ±0.16 


-0.04 ±0.07 

-0.24 ±0.11 

0.11 ±0.11 


0.83 ±0.08 
0.79 ±0.12 
0.84 ±0.17 


-0.07 ±0.08 
-0.35 ±0.11 
0.15 ±0.11 



Table 2. Estimates of the rotational signature and anisotropy 
of the Milky Way stellar halo with BHB stars from SDSS (DR9) 
selected according to two different cuts in the (11-3)0 colour (the 
other filters stay as in ll3p : the tighter cut in {u — g)o aims at 
excluding a suspected blue hook. 



lAihara et all l201ll l3: iMunn et al] |2004| ) and all the above 
statements are confirmed well within one sigma. For those 
previous data releases the selection via photometry was 
obtained with a more generous filter in surface gravity: 
logg € [2.0,3.5], accounting for the higher surface gravity 
estimates fr om the SSPP by on average ~ 0.2 dex in DR8 
w.r.t. DR9 JAhnet al.ll20ia ). The care for the high surface 
gravity stars is justified by the risk of contamination from 
other A-ty po stars when \ogg > 3.2 (e.g. blue stragglers, 
ISirko et al, 2004 ) . All the same, applying the more generous 
filter logg € [2.0,3.5] also to the FS12 DR9 sample does 
not alter our results: the rotation measured by the model 
shifts by ~ 5% and the estimators by 1 — 5kms~^, but not 
in a systematic way. With this cut the sample is enlarged 
by ~ 200 - 400 stars for (u - g)o > 1.15 and (u - g)o > 0.8 
respectively. 



4- 1-2 Systematics in the analysis 

We repeat the same analysis applying the more generous 
cut in u — g colour discussed in H3.2I to include what we 
identified with the onset of the blue hook in the colour-colour 
plane. For our sample we find that the effect is negligible: 
the discrepancy in the rotational signature is less than 7% 
(respectively At; < 0.13) and about 4 — 11% (respectively 
A/3 ~ 0.04 — 0.11) in the anisotropy. Results are reported in 
Table [H 

With regard to a possible error in the BHB absolute 
magnitude (see §refsec:absmag), we re-analyse our sample 
accounting for shifts in Mg of order of ±0.20 and find that 
both the rotational signature and the anisotropy are practi- 
cally unaffected: we measure no systematic variations in ei- 
ther ri or P, with the maximum discrepancy being less than 
4% for the rotation parameter and 5% for the anisotropy. 

We then contaminate our sample with 10% main- 
sequence A stars: we measure no appreciable shift in either 
the rotational signature or the anisotropy previously esti- 
mated. 

Given the on-going debate on the v alue of Rp and o f the 
Sun's peculiar velocity (see discussion in lSchonrichl2012r ). we 
test our results by considering the following variations in the 
radial position of the Sun: Ro = 8.0, 8.2, 8.5 kpc and find no 
difference in either the rotation or anisotropy retrieved. The 
same holds for alterations in the U-component of the Sun's 
velocity that span values from 10 km/s up to 16 km/s. What 
does affect the rotational signature estimated through our 
analysis is the Sun's azimuthal velocity: raising it from the 
adopted Vcq = 220 km/s to Uc© = 250 km/s, makes the ro- 
tation signature 77 rise by 0.5 towards pro-rotation. However, 
the significance of this shift is more intuitive looking at the 
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direct estimators: the +30kms~ change in Vcq translates 
into a shift in the rotation estimate of < lOkms"^ towards 
pro-rotation. 



5 COMPARISON WITH DEASON ET AL. 2011 

Dll presented a study on the kinematics of the Milky Way 
stellar halo using SDSS BHB stars from DR7 and retrieved 
a metal-poor counter-rotating halo and a metal-rich pro- 
rotating halo. They identified BHB halo tracers by cutting 
in colour, surface gravity, effective temperature and position 
to remove the disk and the Sagittarius stream as we do in 
H13p . except they also include stars witn^^l: 




3.5 ^ logy < 4.0 



5.1 Discussion of Dll's selection criteria 

The inclusion by Dll of stars with \ogg £ [3.5,4.0] and 
g ^ 18 exp oses them to the r isk of contamina tion by non- 
BHB stars (|Sirko et al.ll2004 iLee et al.1 12008^1. Therefore, 
the rotational signature retrieved with their sample is more 
prone to be biased: according to the literature value they 
buy their larger sample size with a contamination by mis- 
selected stars of ~ 25% (compared to 10% for our sample, 
see section 3) and the inclusion of faint stars makes their 
sample span distances from the GC up to 90 kpc, while we 
cover a sphere of radius 30 — 40 kpc. Therefore the geomet- 
rical biases outlined in §3.1 will be more significant for the 
sample of Dll than for ours. 

We confirm that if we drop our tighter cuts in surface 
gravity and apparent magnitude, we retrieve the results of 
DIFt Also, we try the heuristic exercise of perturbing the 
distances in the sample of Dll by 20% and verify that the 
inclusion of faint stars makes the rotational signature highly 
susceptible to distance errors. 

We assume a misclassification of AAfg = ±0.40, cor- 
responding to a distance underestimate of ±20%^H Fig. [7] 



^^ Note that given the higher average surfac e gravity value i n 
SDSS data releases prior to DR9 by ~ 0.2 dex l lAhn et al.ll2012l) , 
for DR7 the surface gravity cut in II13I I. reads: 2 < log 3 < 3.5 
(see i|4.1.1l for a discussion). 
^■^ With the sample drawn from SDSS DR7, filtered as pre- 



scribed by Dll and assuming Mg 



0.7, wc find the same two 



populations (a metal-poor one with —3 < [Fc/H] < —2 and 
a metal-rich population identified by — 2 < [Fc/H] < 0) and 
estimate /3 and r] for them separately, according to the like- 
lihood ^. We confirm the trend presented by Dll of a pro- 
rotating metal-rich population versus a counter-rotating metal- 
poor one and match their numerical estimates of rj and /3 with 
the exception for the anisotropy parameter for metal-rich stars. 
Dll find iv,l3Uch = (0.8 ± 0.1,0.1 ± 0.1) and (»?, /3)poor = 
(1.4±0.1,0.2±0.1) while we find (??, ^)rich = (0.7±0.1, 0.5±0.2) 
and (ry, /3)poor = (1-5 ± 0.1, 0.3 ±0.1). We suspect this difference 
traces back to their undescribed handling of the singularity at 
L = by Dll and to the 15% difference in sample size. 
15 "\Yg remark that regardless of the theoretical nature of this ex- 
ercise, these values for AMg are not at all unrealistic (see Section 
2). 



Figure 7. Probability density of rj marginalized over /9 for BHB 
stars from DR7 (see 4.1): blue lines for metal-poor stars and red 
lines for metal-rich ones. Full lines represent Dll's results, who 
assume that Mg = 0.7 for each star in the sample. Dashed lines 
show how the rotational signatures change if the distances have 
been underestimated by 20%, dotted lines represent the signa- 
tures associated with a distance overestimates of the same order. 



shows that for distance overestimates (dotted lines) the ro- 
tational signatures are noisier and the discrepancy between 
the metal-poor (blue lines) and the metal-rich (red lines) 
sample is exacerbated, while for distance underestimates 
(dashed lines) they shift towards non-rotation and the signal 
reduces. On the other hand, when we challenge our sample 
with a perturbation in the adopted absolute magnitude of 
similar intensity, we measure no scatter greater than one 
sigma in the rotational signatures previously estimated. It 
seems helpful to remark that distance errors do not change 
the sign of rotation, but rather enhance (for distance over- 
estimates) or shrink (for distance underestimates) the mea- 
sured rotational signature: this is a pure geometrical conse- 
quence for measuring rotation via the component of velocity 
that at the limit of infinite distance is perpendicular to the 
direction of rotation. Therefore, if there is any rotational 
signal a distance bias will infiate it or reduce it, but if the 
system is non-rotating, errors in distance will have no effect 
on the estimated streaming motion. Distant stars are par- 
ticularly sensitive to distance systematics for the effect of 
a bias scales with distance (see Fig. [2|. This explains why 
our sample is so little affected by perturbations in the dis- 
tance scale, while the signature found by Dll is not robust 
in these respects. 

In conclusion, the rotation measured by Dll is driven 
by stars with logg £ [3.5, 4.0] and g ^ 18. When we remove 
these stars, with the same analysis we measure no rotation. 
We note that this exclusion does not bias the properties of 
our sample: for example the distributions of vioa and of Tefi 
are unaffected by the tighter filters. The i; < 18 cut very rea- 
sonably shifts the peak of our distance distribution towards 
a shorter mean, but has no effects on the other properties. 
Further, this does not affect our conclusions since the sig- 
nal claimed by Dll is constant with radius. The removal of 
loggr ^ 3.5-stars diminishes mainly objects from the metal- 
rich sub-sample and far less from the metal-poor sub-sample. 
Since the controversial signature resides on the metal-poor 
side of the sample, the signal should be enhanced by the 
altered proportions, but this does not happen. 
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Besides the concerns associated with the classification 
of these stars, if the signature measured with Dll's more 
generous filter were a global property of the sample consid- 
ered, removing random elements from it would widen the 
spread around the mean, but keep the latter the samqij. On 
the contrary the comparison shows a shift in the means and 
a decrease in the formal errors, confirming the suspicion that 
the objects with logg £ [3.5,4.0] and g ^ 18 are non-BHB 
stars. 



5.2 Consistency between model and data 

To assess whether there is any systematic in the analysis of 
Dll due to the more generous filters in apparent magnitude 
and surface gravity, we run the consistency check introduced 
on their sample (see tj2.3.1|) . In the subsequent subsection 
I5.2.2l we will trace back the detected inconsistency to a likely 
disc contamination. 



5.2.1 Spatial dependence of the rotational signal 

We bin {I, 6)-space to have cells with useful statistics accord- 
ing to the {l,b) distribution of BHB stars within SDSS: we 
attempt a regular grid and perform the analysis in each bin. 
We compute the formal error in each bin, and since they 
are all of the same order we consider the regular grid a fair 
choice. There is of course a geometrical effect one needs to 
account for: the sphericity of the reference frame causes no 
rotation at the poles {b = ±90) and maximum rotation on 
the galactic plane. However, given the significantly big bins 
chosen, this produces a smooth fall of the rotational signa- 
ture measured from the plane up to the poles, rather than 
jumps. From one bin to the next the mean rotation velocity 
is therefore expected to be consistent within the error bars. 
The only place where we can have a real physical disconti- 
nuity in the rotational signature is the region which includes 
the GC. 

Fig. [8] shows the bins in (/,fe)-space with an arrow in 
each bin pointing in the direction of the recovered rotation 
for the metal-poor component. The arrow lengths are asso- 
ciated to the strength of the rotational signature, and the 
numbered values of 77 are reported above the arrow. We see 
that the best values for the rotational parameter are con- 
sistent, within the errors, with a non rotating stellar halo 
{rj = 1). Fig. [9] shows the same analysis for the metal-rich 
population: here the situation becomes pathological as there 
are jumps from pro-rotating to counter-rotating motions. 
Compare for example the bin centered on (135, —68) with 
the one on (135, —23) or the bin centered on (225,23) with 
the one on (225,68). 

From this, we conclude that the model ([l]) is unable to 
fit the sample of Dll and the parameters estimated using 
this model cannot be used to achieve any further under- 
standing of the system. Indeed, if the model is a bad rep- 
resentation of the data, the parameters estimated from it, 
will be driven by the noise that arises from the discrepancy 



^® Random is here understood as statistically non significant and 
the latter consideration is a consequence of the assumption that 
there are no spectral contamination in the population. 




Figure 8. Analysis of the rotation parameter rj on metal-poor 
BHB stars from SDSS DR7 as selected in Dll and performed 
locally in {I, 6)-space. Arrows pointing right indicate pro-rotation 
and vice versa for arrows pointing left. Their magnitude is a vi- 
sual indication of the rotation strength while the actual value is 
reported above them together with the associated formal error. 
Bins are coloured according to the density of stars (reported in the 
bottom-right corner of every bin). The colour scale is normalized 
per each population. 
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Figure 9. Same as in Fig. ^ for metal-rich stars. Note the jump 
in signature between the bin centered in (135, —68) and the one 
in (135,-23). Or between the bin centered in (225,23) and the 
one in (225, 68) bin. 



between the model and the true distribution, rather than by 
the real signal in the data. 

A second consistency check is presented in Appendix C, 
where we show that for a distribution function like lU and 
for fixed /3, the best-fitting value for the entire population 
should be a linear combination of the best-fitting values of 
any sub-samples one breaks it into, where the coefficients 
are the fractional sizes of each sub-sample with respect to 
the total size. Our results satisfy this constraint, but the 
ones reported by Dll do not. From their Fig 4. rjrich = 0.8, 
VpooT = 1.4, 77totai = 1.1. For fixed /3 = 0.1, having 32% 
metal poor stars and 60% metal rich stars, one should have 
T/totai ~ 0.9 (if one considers the 3549 stars) or 1.0 (if one 
previously excludes from the total sample stars with "no 
assigned metallicity" , i.e. ~ 300 stars). 

The U-estimator (see 3.4.1) supports the conclusions of 
the two sanity checks just introduced: on the sample of Dll 
we measure an overall very weak counter-rotation actually 
consistent with no-rotation at all: {v^) = 4± 10 km/s. How- 
ever, with this sample, the radial velocity of the Sun is offset 
by a bit more than 1 sigma: Uq — 1.4 ± 6.7 km/s. 
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5.2.2 Likely disc contamination m Dll's sample 

The prograde rotation measured by Dll on their metal-rich 
sample is linked to the spatial inconsistency of their data- 
model: Fig. [5] shows that the systematic prograde rotation 
of metal rich objects occurs at low latitudes and fades away 
at higher latitudes. A priori this would be consistent with 
either a disc contamination, the lower halo being set into 
stronger rotation by friction with the disc or the debris of 
an accretion event. 

In an attempt to break the degeneracy, we vary the 
lower altitude cut. When we vary the z-cut from 4 to 6 kpc 
(cutting about ~ 10% of the sample) in the Dll sample, the 
vios-estimator confirms a statistically significant drop in the 
rotational signal of metal-rich stars from 34.1±7.4kms~^ to 
23.5 ± 8.2kms~^. At a cut of 14 kpc the iiios-estimator has 
lost its significance by dropping to 12.5 ± 16.6 kms~^. Hence 
we know that basically all the prograde signal is coming 
from stars close to the Galactic disc. How reliable is it? To 
answer this, we look again at the 3D rotation estimators and 
compare to the other, cleaner datasets. 

In the cleaner FS12 sample we find significantly slower 
rotation of 22.2 ± 12.0 kms^^ for 1^1 > 4kpc and a statisti- 
cally insignificant signal of 11. 6± 12.6 kms~^ for \z\ > 6 kpc, 
while the Xll sample, which should have the lowest contami- 
nation with blue stragglers/main-sequence stars, has no pro- 
grade rotation at all (see Table [2]). Hence, with increasing 
data quality the signal diminishes. 

On the other side the 3D estimators do not display any 
rotation at all. In the most strongly rotating case of Dll we 
find: 



(tT (l^los)) 
but {V4,)u-est 



(34±7)kms"\ 
(-2±10)kms"\ 



{v^} = (3 ± 6) kms 



where the values are taken for \z\ > 4 kpc. This inconsistency 
could be the consequece of part of a stream in the sample, 
but even then its rotation wouldn't be necessarly prograde. 
Further, the 3D estimator (v^) shows a hint of rotation in 
the near-disc region and interestingly gets into better agree- 
ment with 7r~^(wios) when we shorten all distances, while 
the [/-estimator stays consistently at 0. All this points to 
the simpler and therefore favourable hypothesis of a con- 
tamination with fainter disc objects. 

It is worth remarking that for the sample of Dll, disc con- 
tamination is likely to exceed the geometric limits of the disc 
itself: the a priori comparison with the sample of IXue et al.l 
l|2011l ) fiagged half of the stars as probable non-BHB. If a 
blue straggler/main sequence star is mistaken for a BHB 
member, its distance is significantly overestimated, so disc 
contamination due to these objects can affect the BHB sam- 
ple well above the actual disc limits. 

In summary even on the Dll sample the 3D estima- 
tors refute any reliable evidence for a clean prograde rota- 
tion, while the difference between the estimators points to 
a contamination by disc objects. This does not mean that 
no prograde tendency in the near-disc region of the Galactic 
halo exists, but it i s not detectable at the accuracy of the 
iDeason et al.l (|201ll ) sample or other samples. 



5.3 Rotation gradient 

5.3.1 In metallicity 

The gradient in rotation that Dll measure between samples 
of different metallicity could at first sight be interpreted as 
support of the claim of a double-component stellar halo, 
even though Dll w.r.t. their result advocate the need of 
redefining the Sun's azimuthal velocity to 250kms~^ and 
impute the rotation gradient to contamination from an ac- 
cretion event|^ Besides the fact that with stringent quality 
cuts we do not measure any rotation discrepancy between 
metal-poor and metal-rich stars (see §4), we investigate if 
there is any further ground in the sample of Dll that could 
support the thesis of a double-component stellar halo. We 
have searched systematically for evidence of bimodality in 
other distributions and can report that no bimodality is ev- 
ident in any of the following planes: [Fe/H] versus any of 
'^ ~ 9i g ~ i~, \og{g), g-mag or v^; v^ versus any of g-mag, 
u — g or r; OY r versus u ~ g. 



5.3.2 In radius 



On the other hand. ICarollo et al.l l|2007l . I2OI0I ') and 
iKinman et al.l (|2007l . l2012l 'l suggest the existence of a double 
halo in terms of real space position: an inner halo dominant 
at 5-10 kpc and an outer (counter-rotating) halo dominant 
beyond 20 kpc. We investigated this possibility both with 
our samples and with the selection of Dll, but find no evi- 
dence for such bimodality in either the radial density profile 
or in the radial distribution of the rotational signature (e.g. 
outer halo counter-rotating and inner halo pro-rotating). 
Fig. [in] shows the mean-streaming motion as a function of 
galactocentric radius: we split the sample into 5kpc-wide 
bins, values being shown every 2.5 kpc so that every sec- 
ond value is independent. The top panel refers to the FS12 
sample and the bottom one to the low c^ stars in Xll (see 
discussion in the following section): to determine the signifi- 
cance of the deviations from a zero constant rotation we run 
a p-test on both trends. Assuming a certain null-hypothesis 
(model) is true, the p-value is the probability to measure a 
deviation from it as large as the one observed. The levels of 
significance to reject a null-hypothesis are most commonly 
taken to be 1 — 5%. The p-value for the trends shown in 
Fig. [TD] confirms that the deviation from a constant zero ro- 
tation is statistically insignificant for all estimators on all 



^^ We remark that drawing from tlie SDSS DR7 catalogue and 
applying the same cuts as in Dll, lead to obtain a slightly dif- 
ferent sample size with respect to Dll. With the colour, gravity, 
temperature and height filter we have 5523 stars (Dll, 5525) and 
when masking out Sgr stream and any close star (r > 10 kpc) 
the sample shrinks to 3500 stars (Dll, 3549). Of these, 2525 fall 
in the metal-rich population (Dll, 2125) and 975 in the metal- 
poor population (Dl l, 1135) even though we are u sing the same 
metallicity indicator l iWilhelm. Beers fc Gravlll999l) . 
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Figure 10. Mean streaming motion as function of galactocentric 
radius for the FS12 sample (top panel) and stars in the XI 1 sam- 
ple with Cj < 0.95 (bottom panel). In one bin out of eleven for 
FS12 and three out of twelve for Xll, the U-est returned a solar 
value in disagreement with the literature: the level of systemat- 
ics in those bins is therefore relevant and the rotation estimate 
cannot be trusted there (dark green points). 



samples[ 



In Table [3] we report the actual estimates of the mean 
streaming motion for two spatial components, one of stars 
closer than 15kpc in galactocentric radius and one of stars 
further than that. All estimators show[^ a halo weakly rotat- 
ing if at all and that there is no appreciable signal between 
the components either spatially or che mically sele c ted. N ote 
that the prograde rotation found by iKafle et al.l (|2013l ) at 
small radii in their metal-rich sample of BHB stars drawn 
from the catalogue of IXue et al.l (|201ll '). is based on iiios in- 
formation only: we find no grounds to support it looking at 



^* In our analysis we use 1 minus the cumulative distribution of 
the x^-distribution to estimate the p- value. Note that because our 
bins overlap by half their size, we effectively have half the degrees 
of freedom we would have if the bins were all independent. 
^^ The discrepancy between the numerical estimates reflect the 
different sensitivities of the estimators to the sample geometry. 



v^{kms ^) 


i;^-estim 


U-estim 


Ulog-estim 


r < 15 kpc 








-2 < [Fe/H] < 
-3 < [Fe/H] s: -2 


-3 ±8 

10 ±11 

-21 ± 13 


7 ±19 
9 ±24 
6 ±32 


22 ±13 
31 ±16 

7 ±21 


r > 15 kpc 








-2 < [Fe/H] < 
-3 < [Fe/H] s; -2 


7±7 
12 ±10 
0±11 


-4 ±11 
-6 ±15 
-1±18 


4±11 
14 ±15 
-8 ±17 



Table 3. We test the hypothesis of a double component stel- 
lar halo. With the BHB sample from SDSS DR9, the model- 
independent estimators of Section 2 measure no significant signal 
both between the kinematics of the closer stars as opposed to the 
furthest ones and when we discriminate for their metallicity. 



the more robust 3D estimators (see tj4.1.ip . Even though this 
could be the kinematic trace of a debris from an accreted 
galactic fragment, it is helpf ul to note that a plausible 10- 
15 kms~^ aberration in uios ()Schonrichll2012h would cancel 
any hint of rotation also in the wios estimators, leaving basi- 
cally unaltered the 3D esti mators. Further, the discrepancy 
in the rotational signal that iKafle et al.l l)2013l l see between 
metal-rich and metal-poor stars, is enhanced both by an in- 
appropriate choice of distance scale, and by the kinematic 
bias afHicting high Cy stars in the sample of lXue et al.l (12011! ) 
(see below). 



CRITIQUE OF THE SIGNATURE IN THE 
XUE ET AL. SAMPLE 



As in the case of iDeason et al.l (|201ll ). the data-model pair 
for the Xll sample fails the consistency check introduced in 
§5.2, so we need to investigate the sample in detail. 

In Fig. [11] we see that Xll's sample strongly enhances 
the high temperature-high b.y region previously unexplored 
by the photometrically selected sample. Also, the relation 
between Toff and 6^. clearly changes in this region w.r.t. lower 
temperatures. Separating the two trends by eye (green line 
in Fig. Ill[) the rotation discrepancy between metal-rich and 
metal-poor stars below the boundary disappears (metal-rich 
: ri = 0.76 ± 0.15; metal-poor: r] = 0.94 ± 0.15). On the 
other hand, the stars above the line show a strong counter- 
rotation as well as a signal of A?) ~ 0.66 ± 0.19 (metal-rich: 
77 = 1.11 ± 0.06; metal-poor: 77 = 1.77 ± 0.18). 

The few stars above the line in the photometrically se- 
lected sample show the same trend. Their strong under- 
representation (20% versus 50% in Xll) explains the dis- 
crepancy towards Xll (left panel in Fig. Ilip . 



6.1 Kinematics of high c-, stars 

This dichotomy can be linked to the steepness parameter c-y , 
which is a line indicator characteristic of luminosity type. 
While it does not directly select on kinematics we see a 
strong systematic difference within the metal-poor stars be- 
tween high and low Cy. In Fig. [12] we plot the de-projected 
1.0. s. velocities as a function of the line indicator c^ (red 
squares) and show the least square fit to 7v~^{vioB){c-y) in 
the two regions (green lines). We do not apply any kine- 
matic cut in fitting the 7r~^(?;ios)-trends, but for the plot we 
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Figure 12. Deprojected l.o.s. velocity for metal-poor BHB stars 
spectroscopically selected by XI 1 as function of the line indicator 
c-y. Only stars with •n-~^(t)ios) G [—300, 300] kms"^ are shown (red 
squares), but the least square fit to the 7r~^(Dios)-trends plotted 
here with green lines have been obtained on the full sample. 



Figure 11. BHB stars that appear in both Xll's sample and 
in the sample photometrically selected via l|13|l (left panels) and 
the full Xll's sample (right panels). In the top panels we plot 
fe-y versus Tcff: the green line indicates what by eye appears to 
be a change in trend. Xll's sample clearly enhances the region 
above the line: in the colour-colour plane (bottom panels) that is 
associable with the blue hook. The colour code reflects the value 
of the C:^ Sersic parameter. 



1)0 ( kms -'- ) 


V 


t)0-est 


U-est 


''los-cst 


c~, < 0.95 










all 
rich 
poor 


0.90 ±0.17 
0.93 ±0.16 
0.87 ±0.19 


-13 ±8 
-13 ±11 
-12 ±13 


-20 ± 16 

-27 ±21 

-11 ±24 


10 ±11 

11 ±14 

8 ±17 


cj > 0.95 










all 
rich 
poor 


1.11 ±0.13 
0.91 ±0.10 
1.4 ±0.19 


-2 ±7 
-5 ±8 
2 ±10 


-7 ±13 
-10 ±16 
-2 ±21 


-10 ±9 

8±11 

-36 ± 14 



Table 4. Estimates of the rotational signature of BHB stars spec- 
troscopically selected by Xll. We find that stars with high c-y 
present a pathological correlation between line cuts and kine- 
matics and reproduce the popular result of a metal-rich pro- 
rotating halo and a metal-poor counter-rotating halo. Unbiased 
stars instead confirm our findings of a weakly pro-rotating or non- 
rotating at all halo and the absence of signal between populations 
with different metallicity. 



limit 7r~^(«ios) G [— 300, 300] kms"^ to make the fit differ- 
ence in c-y better visible. Metal-poor stars with Cy ^ 0.95 
display a significant retrograde signal that is absent in the 
low Oy stars. We note that these stars are not more distant 
than stars with c^ < 0.95; indeed the density distributions 
of the two sub-samples show no appreciable differences. The 
metal-rich stars do not seem to suffer from this systematic. 
We bisect our sample in both in Cy and [Fe/H] and 
for each of the four sub-samples determine the cumulative 
distributions in «ios and 7r~^(iiioa). We then determine the 
Kolmogorov-Smirnov statistics to ask whether the former 
distributions can be drawn from a common underlying pop- 
ulation. Fig. ll3l shows Kolmogorov-Smirnov tests on the dis- 
crepancy between the distributions of wios (upper panel) and 



7i"~^(i'ios) (lower panel) for the four sub-samples. We plot 
the difference in the cumulative probability distributions of 
metal-poor and metal-rich stars in the two Cj bins (low c^: 
magenta line; high c-,: purple line) as well as the difference 
between sub-samples of same metallicity, but different c^ 
(metal-rich stars: red line; metal-poor stars: blue line). For 
each test we show the 5% significance threshold with hor- 
izontal lines of the same colour. The discrepancy between 
two distributions is statistically significant if this threshold 
is crossed. 

Indeed, we find that the most significant deviation oc- 
curs between metal-poor and metal-rich stars with c-f ^ 1.0 
(or, almost equivalently, for b^ > 9.5) and between metal- 
poor stars with different c^. By contrast, stars with c-, < 1.0 
present no kinematic signal between sub-samples of differ- 
ent metallicity. In addition the significance of the signal be- 
tween metal-poor stars in different c-y bins has exactly the 
same structure as the one between metal-rich and metal- 
poor stars with high c-y (and same significance). Therefore, 
the Kolmogorov-Smirnov test confirms that the signature in 
the Xll sample is indeed caused by high c-y metal-poor stars. 
This inconsistency within the metal-poor sub-sample, raises 
severe doubts about the reality of the observed differencesjfj 

Even though with the current data we cannot rule out 
the possibility that these stars are associated with a real sub- 
structure in the halo, the uios-estimator (see Section 2.2.1) 
measures a non-zero gradient on these stars, pointing again 
to a pipeline problem. 

These stars (~ 500) are hot and very metal-poor so 
their spectra harbour nearly no metal lines that can sup- 
port a radial velocity measurement, so this anomaly occurs 
in exactly the group of stars where we have the least trust 
in the velocity measurement. As can be seen from Table UJ 
neither of the two 3D estimators show an appreciable ret- 
rograde signal. Proper motions seem not to contribute to 
the rotational signal and the weak one measured appears to 

^^ When we opt for the pipeline adopted metallicity instead of 
the WBG metallicity, the discrepancy remains and in fact is ex- 
acerbated. 
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Figure 13. Differences in the cumulative distributions of l.o.s. ve- 
locity (upper plot) and de-projected vio^ (bottom plot) between 
sub-samples of different metallicity in two c-y bins. Red lines rep- 
resent the discrepancy between metal-rich stars with different c-y , 
blue lines the same but for metal-poor stars. The magenta lines 
depict the difference between the velocity distributions of metal- 
rich and metal-poor stars with high Cj , the purple lines the same 
but for stars with low Cy. The horizontal lines are associated 
with the corresponding 5% significance levels of the Kolmogorov- 
Smirnov test. 



come from uios alone. In addition the V-W estimator mea- 
sures an underestimate in distance of / = —0.10 ± 0.09 for 
these stars. Since there cannot be a significant contamina- 
tion with more luminous stars (a contamination with blue 
stragglers or blue hook stars would give the opposite sign), 
this points to unusual kinematics in the vertical direction 
rather than a simple retrograde motion. 

In conclusion, we find that metal-poor stars with 



high 



do show a retrograde signal and we identify 



them as the source of t he line-of-sight ve lo city distribu- 
tion a n omaly observed bylBeers et al.l (|2012l ): lHattori et al.l 
l|2013h : lKafle et al.l (|2013l ). By their inconsistency with their 
low Cj counterparts we flag this either as a pipeline problem 
or a special structure rather than as a second smooth compo- 
nent |fj If it were real, this structure would predominantly 



We also note t hat t he discrepancy in metallicity distribution 
that I Beers et al.l 1 120121 ) advocate as evidence for a double-halo 
arises when comparing stars closer and further than lOkpc. The 
shift of 0.3 dex is purely due to stars between 5 and 10 kpc, which 



contain only hot BHB stars and have large systematic W 
velocities across the sky while not showing a consistent ret- 
rograde motion. On the other hand, trustable stars (low Cy) 
confirm our previous conclusions of a weakly pro-rotating 
or non-rotating halo irrespective of metallicity. Equivalent 
behaviours emerge when splitting the sample according to a 
line cut associated with fo-v. 



7 CONCLUSIONS 

We have analysed the kinematics of BHB stars drawn from 
SDSS/SEGUE (DR7, DR8, DR9), for which the only rehable 
measured component of velocity is uios , using both a selec- 
tion via photometry and the sample of IXue et al.l l|201ll ). 
Using a very simple model and three model-independent es- 
timators we find no evidence for rotation. Also, we measure 
no discrepancy in ro tation between metal-poo r and metal- 
rich stars. Following ICaroUo et a!] l|2007l . |2010| ) we look for 
a discrepancy in the kinematics of the inner and an outer 
halo, but find none in any of the possible selections. In both 
samples we trace back previously claimed rotational differ- 
ences between different metallicity bins to stars that don't 
pass stricter quality criteria. 

Our analysis is not restricted to the sim ple double 
power-law model used by iDeason et al.l (|201ll ). but relies 
on three direct model-independent estimators: one uses the 
de-projected vios and two make use of the full 3D motion. 
The latter require knowledge of proper-motion and while 
individual BHB stars are too remote to have sufficiently ac- 
curate proper-motion estimates, the errors being compara- 
ble to the signal, on the entire sample proper motions prove 
useful by large-number statistics. We checked that the 3D- 
estimators are robust both against change in the used Sloan 
data release and against systematics of t he magnitude of 
the correction derived in ISchonrichI (|2012r ). Systematics of 
higher magnitude would bias the SD-estimators differently: 
the fact that they agree and that we recover the correct 
radial reflex motion Uq confirms that the systematic floor 
is within the formal errors and the proper- mot ions useful. 
Moreover, the use of the 3D-estimators is crucial to test the 
results from the ID-estimates that are highly vulnerable to 
pipeline systematics. A systematic shift in line-of-sight ve- 
locities causes a larger shift in the measured rotation, since 
the errors are amplified by the geometric factors and only 
partly balanced by having both sides of the sky. Further, a 
non-zero "rotation" in the one-dimensional estimators does 
not prove a true retrograde or prograde net motion, but can 
be caused by failure of the assumption of zero mean U and 
W component. 



still have a non negligib le probability of bei ng disk contamina- 
tion (hence why both in lDeason et al.l 1 12011^ and here they are 
tossed out). The ta ngentially biased metal-poor halo found by 
iHattori et al.l 1120131) on the other hand, is a consequence of the 
pseudo-rotation of high c-y stars, but the robustness of their mea- 
surement is also undermined by neglected treatment of possible 
contaminations. A minuscule fluctuation arising from a distance 
misestimate for example, can be enhanced to arbitrary level pro- 
foundly biasing the inferred anisotropy, hence the duty to quantify 
their vulnerability to distance errors and to deviations from the 
assumption of sphericity. 
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Model-wise, overestimates in distance cause velocity 
space to shrink so certain stars are pushed out beyond the 
escape velocity with the consequence that the correct model 
assigns them zero probability. In this case the stars are elim- 
inated from the group that contributes to estimates of the 
population's kinematics. This effect is not mirrored for dis- 
tance underestimates, making the bias more severe. Unless 
observational uncertainties are appropriately folded with the 
model, again the 3D estimators are a necessary confirmation 
of any claim of net rotation. 

We have proposed two general and one model-specific 
a posteriori consistency tests. The main idea being that a 
true rotational signature in the sample should be seen con- 
sistently in any chosen portion of the sky. Re-analysing the 
photometrically selected sample of Deason et al. (2011) we 
can recover their difference between prograde metal-rich and 
retrograde metal-poor populations. However, the signal is in- 
consistent across the sky and among the different estimators. 
In particular, the weak prograde signal in the metal-rich sub- 
sample shrinks considerably when we increase the geometric 
disk cut from \z\ > 4 to \z\ > 6 kpc. From the point of the- 
ory, it is plausble that due to accretion or frictional effects 
with the disc a confined region in the Galactic halo could 
display a mildly prograde rotation. However, the observa- 
tions point to a disc contamination, given that the signal 
strongly weakens for the cleaner samples and that the 3D 
estimators do not confirm the rotation. We note that the 
additional constraints only reduce the level of contamina- 
tion, while they do not a priori bias the inferred kinematics; 
when we repeat the same analysis on XI 1 we measure in 
fact homologous trends (see below). We further argue that 
the anisotropy measured via the adopted model is unreliable 
due to the unphysical L~^^ term in the distribution func- 
tion, which for /3 < produces a bimodal distribution in the 
azimuthal velocity component. 

We r epeat our analysis on the spectroscopically selected 
sample of IXue et al.l (|2011h . On the intersection of the two 
samples we do not find any hint for rotation. However, on 
the entire XI 1 sample we detect a weak offset in v^ be- 
tween metal-poor and metal-rich stars. While there is no 
detectable difference between the velocity distributions of 
metal-poor and metal-rich stars at low c^, we identify the 
metal-poor stars with high c^ as the source of this discrep- 
ancy. A Kolmogorov-Smirnov test reveals that the uios dis- 
tribution of these stars is not only offset significantly from 
that of metal-rich stars, but also from that of their metal- 
poor counterparts with low c^. Furthermore the retrograde 
behaviour in vios is not confirmed by the 3D estimators, and 
anomalous W velocities are detected. In fact these system- 
atics are confined to metal-poor hot BHB stars that have 
very weak to non-detectable metal lines, which are crucial 
for reliable t)ios determination. 

There are two possibilities: either this is a pipeline prob- 
lem or these stars belong to a very metal-poor population 
containing only hot BHB stars. A systematic shift in line-of- 
sight velocities is not unlikely, given that the determination 
of line-of-sight velocities is very difficult in the smooth and 
nearly metal-line free spectra of hot, metal-poor BHB stars 
and that the average motion of stars towards the Galactic 
North and South pole indicates a significant pipeline system- 
atic that increases towards blue, metal-poor stars in general. 
If the component were real, it would not be in smooth retro- 



grade motion, but rather likely form a dispersed stream-like 
structure with large systematic W velocities. High resolution 
spectral analysis of the suspected culprits will be required to 
decide if the SEGUE tJios are trustable for those stars, and 
if so, to identify the peculiar abundance pattern that can be 
expected for such a unique metal-poor structure. 

We tested our results against a variety of sources of un- 
certainty such as selection criteria, absolute magnitude es- 
timate, contamination by main-sequence A stars^the Sun's 
parameters, the effect of reddening from the dislo and dis- 
tance errors up to 20%. The results hold and prove to be 
susceptible only to variations in the Sun's circular veloc- 
ity, for which we provide none the less a trend: changing 
VcQ — 220 km/s to Vc© = 250 km/s shifts the rotational 
signature towards pro-rotation by less than 10 kms~^. 

The failure of the model used in this work teaches us 
that significant efforts have to be devoted to the creation 
and selection of appropriate models and strategies for fit- 
ting them to data. So far in fact, there is no a priori way 
to guess which characteristics a model should possess to fit 
satisfactorily the data and in fact, there is no agreed general 
definition of goodness of fit. There are ways to compare a 
model's ability to represent a certain sample (e.g. by com- 
puting the likelihood of the data given the model), but even 
with the best model among the ones available, there is no 
way to assess its intrinsic uncertainty. 

We believe this area is still largely unexplored and 
should be considered the key step to cover for any mean- 
ingful analysis to be pursued. 



ACKNOWLEDGMENTS 

We thank James Binney and Heather Morrison for useful 
discussions and helpful comments on a draft of this paper. 
We further thank Alls J. Deason and Wyn N. Evans for 
providing the details of their analysis in Dll. F.F. acknowl- 
edges financial support from the Science and Technology 
Facility Council (UK) and from Merton College, Oxford. 
R.S. acknowledges financial support by NASA through Hub- 
ble Fellowship grant HF — 51291.01 awarded by the Space 
Telescope Science Institute, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., for 
NASA, under contract NAS 5-26555. Funding for SDSS-IH 
has been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Science Foundation, 
and the U.S. Department of Energy Office of Science. The 
SDSS-IH web site is http://www.sdss3.org/ 



REFERENCES 

Ahn C. P., et al., 2012, ApJS, 203, 21 
Abazajan K. et al., 2009, ApJS, 182, 543 
Adelman-McCarthy J.K. et al., 2006, ApJS, 172, 634 
Aihara H. et al., 2011, ApJS, 193, 29 
Aihara H. et al., 2011erratum, ApJS, 195, 26 

^^ We studied the effect that the choice to cut for |2| < 4 kpc 
has on our results. We consider different cuts at \z\ = 1, 2, 3, 5, 10 
kpc and observe no change in the signatures retrieved with our 
sample. 



© 2013 RAS, MNRAS OOO.[lH20l 



18 Fermani & Schonrich 



2001, MNRAS, 328, 



Hanson R. B. 



Armandrofl T. E., 1989, AJ, 97, 375 
Baev P. v., Markov H., Spassova N., 

944 
Beers T. C, et al., 2012, ApJ, 746, 34 
Beers T. C, Chiba M., Yoshii Y., Platais I. 

Fuchs B., Rossi S., 2000, AJ, 119, 2866 
Bell E. P., et al., 2008, ApJ, 680, 295 
Belokurov V., et al., 2006, ApJ, 642, 137 
Binney J. J., Merrifield M., 1998, Galactic Astronomy 

(Princeton, Princeton University Press) 
Binney J. J., Tremaine S., 2008, Galactic Dynamics (2d 

ed.; Princeton, Princeton University Press) 
Bullock J. S., Johnston K. V., 2005, ApJ 635, 931 
Byrd G. G., Saarinen S., Valtonen M. J., 1986, MNRAS, 

220, 619 
CaroUo D., et al., 2007, Nature, 450, 1020 
CaroUo D. et al, 2010, ApJ, 712, 692 
Chalonge D., Divan L., 1973, A&A, 23, 69 
, 2000, AJ, 119, 2843 
J., Hewett P. C., et al., 2002, MN- 



Chiba M., Beers T. C 
Clewley L., Warren S 

RAS, 337, 87 
Cooper A. P., Cole S 

406, 744 
Deason A. J., Belokurov V 

411, 1480 (Dll) 
Deason A. J., Belokurov V 

416, 2903 
Dong R., Gunn J., Knapp G., Rockosi C, Blanton M 

2011, AJ, 142, 116 
Dorman B., Rood R. T., O'Connell R 

596 

Lynden Bell D., Sandage A 



Frenk C. S. et al., 2010, MNRAS 



Evans N. W., 2011, MNRAS, 



Evans N. W., 2011, MNRAS, 



Eggen O. S 

136, 748 
Fermani F 
Freeman K. 
Frenk C. S., 
Fukugita M 



W., 1993, ApJ, 419, 
1962, ApJ 



R., 



Schonrich R., 2013, MNRAS, 430, 1294 
, 2011, arXiv:1108. 5028 1^1 
, White S. D. M., 1980, MNRAS, 193, 295 
Ichikawa T., Gunn J. E., Doi M., Shimasaku 
K., Schneider D. P., 1996, AJ, 111, 1748 

Hartwick F. D. A., 1987, NATO ASIC Proc. 207, 281 

Hattori K., et al. 2013, ApJ, 763, L17 

Heavens A., 2010, arXiv:0906.0664v3 

Kafle P. R., et al., 2013, MNRAS, 430, 2973 

Kaczmarczik M. C, et al., 2009, AJ, 138, 19 

Kinman T. D., Cacciari C, Bragaglia A., Smart R., 2012, 
Spagna A., MNRAS, 422, 2116 

Kinman T. D., Cacciari C, Bragaglia A., Buzzoni A., 
Spagna A., 2007, MNRAS, 375, 1381 

Kinman T. D., 1995, lAUS, 164, 75 

Klement R. J., 2010, A&ARv, 18, 567 

Layden A., Hanson R., Hawley S. L., Klemola A. R., Han- 
ley C. J., 1996, AJ, 112, 2110 

Lee Y. S., et al., 2008, AJ, 136, 2022 

Lee Y. S., et al., 2008, AJ, 136, 2050 

Majewski S. R., 1992, ApJS, 78, 87 

McMillan P. J., Binney J. J., 2012, MNRAS, 419, 2251 

Munn J. A., et al., 2004, AJ, 127, 3034 

Munn J.A., et al., 2008, AJ, 136, 895 

Murante G., Poglio E., Curir A., Villalobos A., 2010, ApJL, 
716, 115 

Nordstrom B., et al., 2004, A&A, 418, 989 

Norris J. E., Ryan S. G., 1989, ApJL, 336, 17 



Odenkirchen M., Grebel E. K., Rockosi C. M., Dehnen W., 

Ibata R., et al, 2001, ApJL, 548, L165 
Philip A. G. D., Adelman S. J., 1994, Workshop on 

Database for Galactic Structure (Schenectady, p. 245) 
Pietrinferni A., Cassisi S., Salaris M., Castelh F., 2004, 

ApJ, 612, 168 
Pietrinferni A., Cassisi S., Salaris M., Castelh F., 2006, 

ApJ, 642, 797 
Press W. H., Schechter P., 1974, ApJ, 187, 425 
Preston G. W., Shectman S. A., Beers T. C, 1991, ApJ, 

375, 121 
Quinn P. J., Goodman J., 1986, ApJ, 309, 472 
Ryan S. G., 1992, AJ, 104, 1145 
Samland M., Gerhard O. E., 2003 , A&A, 399, 961 
Sandage A., 1970, ApJ, 162, 841 
Schneider D.P. et al., 2010, AJ, 139, 2360 
Schonrich R., Binney J., Asplund M., 2012, MNRAS, 420, 

1281 (SBA12) 
Schonrich R., Binney J., Dehnen W., 2010, MNRAS, 403, 

1829 
Schonrich R., Asplund M., Casagrande L., 2011, MNRAS, 

415, 3807 
Schonrich R., 2012, MNRAS, 427, 274 
Searle L., Rodgers A. W., 1966, ApJ 143, 809 
Searle L., Zinn R., 1978, ApJ 225, 357 
Sesar B. et al., 2010, ApJ, 708, 717 
Sirko E., et al., 2004, AJ, 127, 899 
Smith M. C, et al., 2007, MNRAS, 379, 755 
Steinmetz M., 2003, in GAIA Spectroscopy: Science and 

Technology, ed. U. Munari, ASP Conf. ser., 298 
Steinmetz M., Zwitter T., Siebert A., Watson F. G., Free- 
man K. C, et al., 2006, AJ, 132, 1645 
Trumpler R. J., Weaver H. F., 1962, Statistical Astronomy 

(New York, Dover Publications) 
Verde L., 2007, arXiv:0 712.3028y 2 
White,S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
Wilhelm R., Beers T. C, Gray R. O., 1999, AJ, 117, 2308 
Xue X. X., Rix H. W., Zhao G., et al., 2008, ApJ, 684, 1143 
Xue X. X., Rix H. W., Yanny B., et al., 2011, ApJ, 738, 79 
Yanny B., et al., 2000, ApJ, 540, 825 
Zinn R., 1985, ApJ, 293, 423 



APPENDIX A: MCMC CHAIN CONVERGENCE 
AND TEST OF IMPLEMENTATION 

1 Sampling parameter space 

We now explain how we retrieve {-q, /3) out of the log- 
likelihood @. 

We construct a MCMC chain in parameter space, where 
the probability density is ([3]). When can such a chain can be 
considered to have converged to a realisation of the model 
probabil ity den s ity? W e establish convergence using the cri- 
teria by IVerdel (|2007l '): given T chains of length 2A^, each 
with a different starting point, one estimates the chain-to- 
chain variance B and the average variance of each chain W, 
which is an underestimate of the true variance: 









(1) 
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T(N - 1) ^^ 

J = l i = JV + l 
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A3 



aJ^2 



(2) 



where 9^ is the point in parameter space in position i of 
chain j, 6-' is the mean of the chain j, ^ is the mean of all 
the chains and the burn-in time A^ is half the length of the 
chains. 

An overestimate of the variance is 



V 



2 I 



B 

NT 
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{N ~l)W + B 



'4 



(3) 



where a'^ is the true variance. iHeavend (|2010l ) suggests that 
for convergence one should require -^ ~ 1.01 — 1.1: the preci- 
sion reported is empirical and does not follow from a rigorous 
proof. We will therefore opt for a much more conservative 
threshold. 

Here, to avoid numerical issues when dividing quantities 
that are potentially null, we adopt the following equivalent 
criterion: 



Y^ = {l + e)^V-W = eW, 

where e > and 

1 W 



(4) 



(5) 



T chains are therefore run separately but at the same 
pace: every time a new step is added to the chains, W and 
B are re-computed, and if the convergence criterion above 
is not yet satisfied, the chains keep growing. This avoids the 
need for an a priori fixed chain length and burn- in time. In 
our code we used: T = 2, e = 10~^ and achieved convergence 
for iV typically of order ~ 10^ - 10^. 



2 Test of implementation 

We test the implementation of likelihood analysis just in- 
troduced against 12 mock samples: we generate a sample 
by MCMC sampling phase-space (4 x 10 points) according 
to the distribution function in equation U]) and random pa- 
rameters (7?,/3) G [0,2] X [-0.5,0.5]. The reliability of the 
sample obtained is checked by comparing the density pro- 
file and rotational signature derived from star-counts of the 
synthetic catalogue produced with the theoretical profiles, 
revealing perfect agreement. 

Subsequently, we determine the probability density of 
the model by MCMC sampling (77, /3)-space using the likeli- 
hood of the pseudo-data given the model. The original pa- 
rameters are recovered within one sigma as demonstrated in 
Fig. [T] From the distribution of the deviations we also infer 
that the error estimates are reasonable. 



APPENDIX B: SAMPLE SELECTION 

The following SQL query was used to drawn data from the 
SDSS archive: 




A,, 



Figure 1. Performance of the MCMC fitting technique: the plot 
shows the discrepancy between the actual parameters and the 
fitted ones, weighted by the formal error. The fit is performed 
by MCMC sampling parameter space using the likelihood of the 
pseudo-data given the model. The green circumference represents 
the one a contour, the red one the two cr contour. 



SELECT 

sp.PLATE,sp.FIBERID,sp.specObjID,ph.rerun, 

,sp.TARGETSTRING,sp.mjd,sp.flag,sp.fehadop,sp.loggadop, 

sp.fehadopunc,sp.FEHWBG,sp.TEFFWBC,sp.LOGGWBG, 

ph.psfmagu,ph.extinctionu,ph.psfmagg,ph.extinctiong, 

ph.psfmagr,ph.extinctionr,ph.psfmagi,ph.extinctioni, 

ph.psfmagz,ph.extinctionz,ph.l,ph.b,ph.ra,ph.dec, 

m.pmL,m.pmB,m.pmRa,m.pmDec,m.pmRaErr,m.pmDecErr, 

m.sigRa,m.sigDec,m.nFit,m.dist22,m. match, 

sp.ELODIERVFINAL,sp.ELODIERVFINALERR,sp.SNR, 

sp.TEFFADOP,sp.TEFFADOPUNC 

FROM sppParams sp 

JOIN PhotoObjAU ph on sp.bestobjid = ph.objid 

JOIN ProperMotions m on m.objid = ph.objid 

WHERE sp.sciencePrimary = 1 

AND ph.CLEAN^l 

AND sp.TARGETSTRING = 'BHB' 

The masking of Sgr is performed in right ascension- 
declination space (q, 5) according to the following polygon 
kindly provided by the authors of Dll on request: 

ai = [245,225,211.5,210,208,200.5,199.3,197.5,197, 

180, 160, 140, 110, 113, 140, 160, 180, 202, 225, 245], 
5i = [-3,8,14.2,12.5,16,19.5,16.5,16.9,21.5,28,33, 
35, 38, 15, 13.5, 11, 6, -3, -3, -3]; 
U a2 = [10,50,70,30, 10], ^2 = [-10,20, 10, -20, -10]. 



APPENDIX C 



,X(N)} 
Xm} U 



Theorem: Let the sample of stars be X = {3::{i) 
divided in the sub-samples X = Xi U X^ = {xi, 
{a;^, . . . ,a;"}. If P{'r)\X) is the probability of the model rel- 
ative to the parameter 77 given the data and ritest,x is the 
value that maximize it for the type of distribution function 
under consideration: 
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card{Xi) card{X^) 

Vbest,X — ,,^^ rjbest,Xi + ,,^^ Vbest.X'^ 



card{X) 



card{X) 



Proof: 

Since / oc 1 + (1 — 77) tanh(-^), assuming A <^ L^ and at 
fixed /3 (such that the normahzation constant, independent 
of ry, is Uc): 

Piv\X) =\\P{v\x(i)) = — ^ '-^ 

Tic 
i = l 

where H{x) = 2H{x) — 1 and H{x) is the Heaviside 
function. Then one can split the above in the positive and 
negative contributions, obtaining: 



^(v\x)^l[l-R- 



where N^ is the number of stars with L^ > and the 
opposite for N~ . The probabiUty is maximum for — pp — - = 

i^ . The same holds for the 



0, which happen for 1 — r; = ■^— ^ 
two sub-samples of X: calling ni"*" the number of stars with 
positive Lz that belong to Xi and n"*" the ones belonging to 
X and similarly for stars with negative L^, we have: 



Vbest,Xi — 



2m" 



and 



Vbest,X^ — 



2n" 



m n 

with A^"*" = m'^ +n'^ and A'^" = m~ +n~ . We conclude, 
showing that rjbeat.x can be expressed as follows: 



Vbest.X 



2(m + n 
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m 2n 

N'^~^ ' N 



This completes the proof. 
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